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1 INTRODUCTION 



In two recent papers Elias et. al. showed that by analyzing the sequence 
of occurrence of the first zeros of the denominator and numerator of a Pade 
approximant beta function, one can separate out two possible behaviors of 
Pade QCD couplant in the infra-red (IR) region. The couplant may exhibit 
an IR attractor behavior of Kogan-Shifman type |^ in which the PQCD cou- 
plant bifurcates at some critical cut-off momentum into a second (upper) 
branch, in a manner that leaves the infra-red region < n < Hc totally free of 
any PQCD color force and its bifurcated branch. We may call this scenario 
II. The alternative scenario I is that the PQCD couplant freezes at the critical 
momentum point /i = fic, being its stable infra-red fixed point, and remains 
so frozen for all < /i < /^c in the manner of Mattingly-Stevenson Q . 

The finding of Elias et. al. by their denominator- numerator zero ap- 
proach is that scenario II holds for all < Nj < 5, regardless of which Pade 
approximants, [1|2], [2|1], [1|3], [3|1], or [2|2] they used. On the other hand, 
they found that scenario I holds for all 9 < Nf < 16. For the flavor numbers, 
Nf = 6, 7, 8, their method led to indeterminate results. Their positive de- 
nominator zeros ceased to exist for Nf > 6, but then their numerator zeros 
were found to be negative for A^j = 6, 7, 8. As a result their method did not 
resolve the above question for these flavor states of Pade QCD. 

We have re-examined the above question, using a graphical method to 
analyze the Pade QCD couplant equation of evolution, and applying this not 
only to the NNNLO [2|1], [1|2] Pade approximants studied by Elias et. al. 
but also to the optimized [1|1] Pade approximant more directly related to 
the optimized NNLO truncated PQCD beta function studied by Mattingly 
and Stevenson. 

While confirming some of the findings of Elias et. al., we found sub- 
stantially new features. In particular, we found for all flavors, not just the 
Elias et. al. bifurcation at momentum point Q = fJ^ = fJ^c = Qmm, but a 
second bifurcation at Q = fi = Qmax, where in all cases Qmax > Qmin- More 
explicitly, we found the Pade couplant equation yielding three independent 
solutions or component couplants ai, 02, 03 that are however joined together 
into a continuous chain-like spiral structure at the two momentum points 
Qmin and Qmax that can be called couplant bifurcation points. 

The 03 solution arising from the second bifurcation, is found to increase 
towards lower energies, running past the original PQCD bifurcation point at 
Q = He, into the region < fic, the implication being that this second 
bifurcation upper branch is either a genuine NPQCD component, or else a 
persisting Landau pole structure coming from the original truncated PQCD 
beta function on which the Pade beta function was based. 
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By noting from our graphical structures and plots, that the momentum 
gap < Q < Qmin is not zero for all 9 < Nf < 16, we found that scenario II 
actually holds also for the flavor states 9 < Nf < 16 as for all < Nf < 8, in 
difference to the findings of Elias et. al. The only amplification is that one has 
to go to sufficiently low momentum to see the persisting infra-red attractor 
points at /i = /ic = Qmm- This apart, the gap < /i < /ic remains finite and 
empty of any PQCD couplant even for the highest flavors Nf = 16, 15, 14... 
The result is our new feature that there is intrinsically in Fade QCD, no 
couplant freezing for all flavors < Nf < 16. 

We present the above analysis and findings as follows. In section 2, we 
collect together the basic QCD couplant equations we need, together with 
their Fade approximants, and the optimized formulation of the particular 
NNLO Fade beta function. In section 3 , we present our numerical and 
graphical method of solving the Fade QCD couplant equations of orders 
[1|1], [2|1], [1|2] that we considered. Then in sections 4 and 5, we state the 
main features and results we found with these Fade QCDs. In section 6, 
we use our graphical method to show explicitly that infra-red scenario II 
holds intrinsically for all flavors < Nf < 16. In section 7, we relate our 
results and findings to the IR fixed point results of Banks and Zaks , and 
Stevenson et. al. [^, ^ ^, |[. Our summary and conclusions are stated in 
section 8. 



2 BASIC QCD COUPLANT EQUATIONS 
AND THEIR FADE APPROXIMANTS 

2.1 The basic dynamical Equation of evolution of QCD 
coupling Constant 

The starting point is the basic RG dynamical equation of evolution of QCD 
coupling constant given by: 

/^^^ = (1) 

where /i is some arbitrary energy renormalization scale, and is 
the (abstract) QCD beta function. Eqn. (|I]) as it stands is believed to have 
a universal validity for all of QCD, that is to hold over a wide range of 
values of momentum scale /i, except if the beta function is later specially 
parameterized or truncated. In the asymptotically free (AF) region of QCD, 
the beta function is parameterized in the perturbation series form: 
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(3{a) = — 6a^(l + ca + 020^ + c^a^ + 040^ + ■ ■ ■ ^ 00) (2) 
where a is the QCD couphng constant in the form a = a^/Tr, while the 



first few expansion coefficients 6, c, 02,03 have the specific values P, 10, 11 



33 - 2Nf 
6 



2/9o (3) 



153- 19Ar^ ^ ,^ 
2(33 -2^ =^^/^° 



C2{MS) 



16(33 - 2Nf 



2857 5033 _ 325 



Nf + —N^f 



18 ' 54 



(5) 



as well as the four-loop term C3{MS) computed not so long ago by Rit- 



bergen et. al [T^, and given by C3{MS) = (^3/ (3o, where 



/33 = 114.23033 - 27.133944A^/ + 1.5823791A^; + 5.85669582 x 10~^iV; (6) 

The rest of the coefficients q, i = 4, 5, 6, . . . cxd are unknown so that QCD 
in making any direct use of eqn (|^) can at most use only the following trun- 
cated beta functions; 

1. The leading order form (LO), where P{as) is truncated to: 

/3(a) ^ -6a2 = /3«(a) (7) 

2. The next-to-leading order form (NLO), where (3{as) is truncated to: 

(3{a) ^ -ba\l + ca) = /J^^) (a) (8) 

3. The third order form or the NNLO truncation where we approximate 
the QCD beta function by : 

p{a) ^ -fea^ (1 + ca + Caa^) = p^^'^ (a) (9) 

4. The fourth order form or the NNNLO truncation where we exhaust all 
presently known terms of eqn. (^ and approximate l3{as) as: 

P{a) ^ -ba'^il + ca + c^a^ + csa^) = /?(^)(a) (10) 
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These truncations now invariably restrict the range of vahdity of eqn. (|I|) 
to only the purely PQCD region. 

The Pade approximant for QCD is a means of extending any one of these 
truncated series into a specific infinite series of Fade form, that is hoped to 
correct to some extent, the effects of truncation on QCD dynamics given 
by eqn. ([|), compared to any direct use of the pure truncated eqns. (|^) 
to (|l^). Such Fade extensions are then substituted into the general solution 
of eqn. (|1]), given in one form as: 



T = blnS- 
A 



dx — b 



°o dx 



a P(^){X) 



where A is related to Aqcd by 

Aqcd 



In 



A 



c, 2c 



If we substituted eqns. (|^ to (H) in succession into eqn. (IT] 
obtain FQCD couplant equations of the form: 

r = 6 In S = — 
A a 



T = 


61n ^ = 
A 


1 

- ^ 

a 


Q _ 
A 


- + - In 

a 2 


\X\ 



+ cln 



ca 



c 



1 + ca 
2 - 2c, 



2co 



2ac2 + c 



(11) 

(12) 
we will 

(13) 
(14) 
(15) 



where X = C2a^ + ca + 1. 

The main feature of these couplant equations is the existence of one cou- 
plant solution at any given high momentum Q, (or /i), with this solution 
going to infinity when the momentum decreases to a critical value given by 
Q = Aqcd, analogous to a Landau pole. The behavior cuts off any access 
from the FQCD region, into the infra-red region < Q < Aqcd of QCD. We 
are now to examine how this picture of FQCD IR region changes when we re- 
place the above truncated FQCD beta functions by their Fade approximants 
and couplant equations. 



2.2 The Pade approximant QCD. 

According to the Fade principle to Wjj, given a generic infinite series: 



(16) 



n=0 
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whose terms are unknown except the first few p terms, meaning the given 
series is in effect a truncated series of only p usable terms, one can construct 
an infinite series approximation for the same S(x) but now of fully known 
terms, by solving the following identity equation : 



S{N\M){x) = J2 '^nX"' (17) 



n=0 



where explicitly 



S{N\M){x) 



1 + liix + 7i2X^ H h 'Jinx 



N 



1 + l2ix + 7223^^ H h 72Ma;*^ 

with N and M being integers chosen such that N + M = p. The quantities 
jij are unknown coefficients which can all be determined in terms of the p 
known coefficients Cn of eqn. (|1^). 

The quantity S(^n\m){x) as explicitly given in eqn. ([18D is what is called 
the Pade infinite series approximation (approximant) of the original infinite 
series S{x) given in eqn. (|l^). By inverting the denominator and expanding 
out with the numerator, the Pade function is seen to be an infinite series 
of fully known terms. By construction, its first p terms are identical in 
value with the first p terms of S{x), but beyond these first p terms, other 
corresponding terms in general differ in their values were one to compute the 
unknown terms of S{x) and compare with the fully known terms of S(^j\f\M){x) 
as done by Ellis et.al. |]TB[. One finds in general good agreement in the various 
studies to [^] carried out with the Pade approximant. 



Depending on our choice of values of N and M for a given p, the Pade series 
is said to be of order (N, M). Based directly on the only known terms of QCD 
beta function given in eqns. (|^) to @, the only non-trivial Pade approximants 
we can consider, are the (1,1), (0,2) at NNLO level, and the (1,2) (2,1) and 
(0,3) at NNNLO level of Pade summation. In this paper we have analyzed 
the (1,1), (1,2) and (2,1) Pade approximants as representative enough to 
show us the general features of Pade QCD couplant equations. For the (1,1) 
NNLO Pade approximant case, we have also analyzed its optimized form, 
providing further information on the nature of the Pade couplant equation 
of evolution. We give details below of these Pade approximants and couplant 
equations. 

1. The [1|1] Pade approximant and couplant equation 

At the NNLO order the [1|1] Pade approximant is given by: 

•^(i|i) « = T— 19 

1 + 72a 
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with the corresponding Pade beta function being : 

/?(3^')(a) = -ba'^Snmia) = -ba' ( 

\ 1 + 72a 



(20) 



where 71 and 72 are the two unknown Pade coefficients to be determined 
from the identity equation: 



ba (1 + 071) (1 + 072) = —ba (1 + ca + ) + 



00 (21) 



We deduce that: 



C2 

7i = c 

c 



72 = — 



C2 



Then eqn. (P) becomes replaced by: 



(22) 
(23) 



/3 



(3)/ 



/3 



(3P)/ 



2 / 1 + 071 



-ba 

\ 1 + 072 / 

-6a^(l + a7i)(l + a72)"^ 



(24) 



Substituting eqn ( pOD into eqn. (pAj), we obtain the dynamical relation- 
ship between the Pade QCD couplant a, and the dynamical variables 
Q and Nf of QCD. This dynamical equation of evolution of as{Q, Nf) 
becomes after explicit integration: 



T = bin — = — hem 
A a 



a7i 



1 + 071 



(25) 



This (1,1) Pade QCD couplant equation (p5D can be analyzed for its 
features particularly in the infra-red region and compared with the 
Landau pole structure of truncated eqn. (^. 

If we were to analyze this NNLO Pade QCD by the method of Elias et. 
al. 0, 0] , we would focus on the Pade beta function eqn. ( PU] ) and note 
the sequence of occurrence of its denominator and numerator zeros. 
The sequence is shown in Table |I|, and would lead one to expect, a 
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Kogan-Shifman type infra-red attractor behavior for < Nf < 8, 
while the Mattingly-Stevenson |^ infra-red frozen couplant behavior 
may exist for 9 < Nf < 16. We will see later that analyzing eqn. by 
our graphical method described below, provides a lot more information 
on the above scenario, and alters the conclusion significantly. 

The [2|1] Pade approximant and couplant equation 

We write down also the [2|1] Pade approximant given by: 

Pl2\i] W = -ba ^(2|i)(a) = -ba — (26) 

\ J- + ^21«. / 

where 



721 
722 
^21 



C — 



C2 



C3 

C2 
CC3 

C2 



C3 
C2 



(27) 



Substituting this Pade beta function into eqn. ( pJ]) we obtain the [2|1] 
Pade QCD couplant equation: 



r = 61n^ = _+ (^2i-r/2i)ln^ 
A a A I A I 

_^ 121 - 2722 - ^21721 

2722a + 721 
+ r/2i-72i + ^ (28) 

721 

where X = a? 122 + 0721 + 1- 

The [1|2] Pade approximant and couplant equation 

In the same way, the [1|2] Pade approximant is given by: 
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where 



7ii = 




- 2CC2 + C3 










- C2 






- CC2 


Vii = 








C2 


- C2 






- CC3 


rii2 = 








C2 


- C2 



(30) 



Substituting this Pade beta function into eqn. (|TT]) we obtain the r 
equation for the [1|2] Pade QCD: 

T = b\nS- = - + (?7n -7n)ln|a7ii + 1| 
A a 



iVu - 7ii) In |a7ii| 
V12 

7ii 



^ln|a7n + l| (31) 



2.3 Optimization of the NNLO [1|1] Pade Approximant 

Regarding the [1|1] Pade approximant, we can come closer to the work of 
Mattingly-Stevenson for later comparison purposes, by subjecting our [1|1] 
Pade approximant to the PMS optimization principle ^, ^ , and obtaining 
in the process a slightly modified eqn. (^). The details are as follows. 

The PMS principle derives from the fact that a physical observable such 
as the cross section for (e+e^ hadrons), parameterized as an infinite 
perturbative series, has the same physical value irrespective of the RS used to 
compute its perturbative coefficients. However if the series is truncated, (with 
a corresponding truncation of the beta function), the physical observable 
becomes RS dependent. The idea of Stevenson is that one can minimize 
this dependence of a truncated physical observable on its RS variables by 
using only those values of the RS variables that satisfy his equation (called 
optimization equation) : 

diRS) ^^^^ 

where (RS) here stands for the complete set of RS variables on which the 
truncated physical observable R^^-^Q), (and its correspondingly truncated 
beta function P^'^\a)) depend. By solving the set of simultaneous equations 
implied by (|3^), one obtains the optimal or optimized set of values of the RS 
variables one should use in the truncated equations for R^^^ (Q) , together with 
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a corresponding optimized QCD couplant a and its truncated beta function, 
all involved in i?*^") (Q) . In the case of eqn. (^) taken along with some third 
order truncated physical observable R),+^.{Q), the optimization equations we 
need to set up are: 



-ha' 



1 + a7i ' 

1 + ^72, 



(33) 



(34) 



where 71 and 72 are as already given in eqns ( p^ ) , (^). Also we have 
from eqn. (|l]), in the form of eqns. (|ll]), the basic dynamical equation for 
«s(Q): 



r = om— = — hem 
A a 



ail 



1 + a7i 



(35) 



Identifying r and C2 as the only RS variables present in eqns. ( |5BD , ( |5^ ) , (§^), 
the optimization eqn. ( p2D takes the form of two equations: 



'A I 



b da 



i?(=^Hg) = o 



(36) 
(37) 



where /?(2)(a) = ^. 
Substituting eqns. ( 



, (|3^ ) into eqns. (|36| ) , (pT]), we obtain from (pB]): 



(1 + 072) + (1 + 072) (2ari + ac) - (1 + 071) (1 + 2ari + 3aV2) = (38) 
Next, from eq. (|37D we obtain: 

- (1 + 072) + (1 + 2ari + 30^2) = (39) 

Equations ( p5| ) , (pSf ) , (^) become the simultaneous equations to solve 
for the optimized quantities: a = a]T = f]r2 = f2',C2 = C2 ; being the 
optimized values for which our physical reference system R^^\Q) would be 
minimally affected by RS changes. 

To the optimization constraints in eqns. (^5|) , ( |5^ ) , (P5|), one now adds 
the intrinsic constraints on RS variables provided by the general existence 
in PQCD of certain combinations of RS variables that are themselves RS 
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invariants 0, ^ In our case there are two such RS invariants, pi 
and p2, given exphcitly by: 

Pi = T — Ti = f — fi = invariant. (40) 



P2 = ^2 + C2 - (ri + -cf = r2 + C2- (ri + -c)^ = invariant. (41) 

where the RS variables on the right hand side are individually dependent 
on the RS used to evaluate them (except b and c), but pi and p2 have nu- 
merical values not dependent on the RS used to evaluate each combination 
of RS variables shown on the right hand side. Because of this RS invar iance 
property of pi and p2, one can obtain a numerical value for them in any 
convenient RS, and these numerical values of pi and p2 become usable as 
input data into the main optimization eqns. ( pSj ) , ( pH] ) , (p9|). Specifically 
we used the MS computed values of pi and p2 given by P, Although 
Pi and p2 now become mere numbers, eqns. (|^) , (^ , (|^) , (^) , (^Tj), 
taken together, show that these numerical values of pi and p2 still depend 
on the values we assign to momentum Q and flavor number Nf . Denoting 
these numerical quantities that depend for their value only on our choice of 
Q and Nf by pi{Q,Nf) and p2{Nf), the complete set of our optimization 
equations becomes eqns. ( p5D , (|38D , (p9D, together with the numerical value 
constraints provided by pi{Q,Nf) and p2{Nf), for the combinations of RS 
variables shown on the right hand side of eqns. (^2|) , (^3]) below: 

pi{Q,Nf) = T-n (42) 

P2{Q,Nf) = r2 + C2-{n + ^cf (43) 

It is these five simultaneous equations (eqns. (|35| ) , (|38[) , (^) , (^2[) , (^)) 
involving five RS variables (r, a, ri, r2, C2), (but r and a are not independent), 
that we now solve by a process of elimination to obtain one final equation 
for our optimized couplant a = a, in terms of optimized f, as our best value 
(or our minimally RS dependent) solution of eqn. (0) into which the NNLO 
Pade approximant eqn. (^) was already substituted in the course of the 
optimization process. 

The elimination process proceeds as follows: we substitute eqns. (|22|) 
and ( PBD into eqns ( |5BD and (jS^) obtaining: 

3r2 + 2riC + C2 + 3r2(c- -)a = (44) 



11 



and 



We then solve eqns. (|4^) 
r2 = f2- The result is: 



+ 2ari + SaVa = (45) 
51) to obtain the optimized values of ri = fi, 



ri 
r2 







(46) 



Substituting these optimized values ri and r2 into eqns. (^^ , (^31) we 
obtain: 



1 fC2 



f) 



r + - — 



P2(A^/-) = -C2 ^ 

as the optimization constraints r and C2 have to satisfy, where from 
eqns. (|l|) , (|20|), r is already given by eqn. (^). We can solve for optimized 
C2 by regarding eqn. (|48|) as a quadratic equation in C2 written as: 



(47) 
(48) 



4 - (16C')C2 + (40^2 + 







We get that : 



C2 = C2 = 3c^ ± c/D = c(3c ± Vd) = 2c(^^c±^Vd 



(49) 



(50) 



where D = 9c^ — (4p2 + c^) = 8c^ — 4p2- If we now substitute eqn. ( ^ ) 
into eqn. (|47|) we get the optimized r = f given by 



T = Pi(g,iv^ 



/J 



3c ± /D 



(51) 



It is this T we finally substitute into eqn. ( P5| ) to obtain our optimized 
Pade couplant a = d that satisfies our universal eqn ([l|) with eqns. ( pOD , (U^ 
dH) , dH). The result is that: 



- + cln 

a 



1 + a7i 



Pi(Q,iV. 



/J 



3c ± /d' 



(52) 



is our equation for determining d of our [1|1] Pade QCD. Next noting 
from eqns. (0) and (^3]) that 71 = c — ^, meaning 71 = c — ^ we can 
rewrite: 
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l + dji = ca-a(3c±/D) + l = 1 + ca - d{3c ± Vd) = l+cd-2dP (53) 
where P = |(3c ± VD), giving : 



- + cln 

a 



ka 



1 + cd- 2dP 



(54) 



where k = c — 2P. 

The final result is that our optimized [1|1] Fade QCD couplant equation 
is given by: 



Pi{Q,Nf)-P = - + c\n 



ka 



1 + cd- 2dP 



(55) 



To simplify notation, we shall from here on, drop the bar over the opti- 
mized couplant in eqn. (|55|), and simply write a in the remaining parts of 
this paper. 



2.4 Pade beta function denominator and numerator 
zero arguments applied to the optimization equa- 
tions 

As in the case of eqn. (^) discussed earlier, before using eqn. ( |55| ) to investi- 
gate the behavior of Pade QCD couplant with momentum scale Q, especially 
near the infra-red region, one can use Elias et. al type argument H, to 
indicate what to expect from eqn. (0). In this case, we apply the Pade beta 
function denominator and numerator zero tests to the complete set of the 
optimization equations given as eqns. (|35| ) , (^) , ( |39D , (|42D , (|i3|). The 
question is do the equations remain consistent and viable if we solve them 
under the specific condition: (1 -|- 072) = with (1 + 071) 7^ 0, which will 
yield infra-red couplant behavior of the Kogan-Shifman type. Correspond- 
ingly we ask if the same equations remain consistent and viable if we require 
to solve them under the condition that (1 + 071) = 0, but (1 + 072) 7^ which 
will yield Mattingly-Stevenson type infra-red couplant behavior. We find by 
looking at eqns. 



(|38|) , (|3|) , dg) , (H), that while these equations 
remain viable under the denominator zero condition, they lose this viability 
under the numerator zero condition. In particular eqn. ( p5| ) blows up for 
(1 + 071) = 0, 

The indication would then be that while eqn. (|S3D is likely to give the 
Kogan-Shifman type infra-red behavior in most cases, the opposite Mattingly- 
Stevenson type behavior will probably not be found for any fiavor state of 
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eqn. (pq). However, totally independent of the above Elias et. al type con- 
siderations, we now present our own graphical and numeraical method of 
analyzing the Pade couplant eqns. , ( P5| ) , ( PT| ) and the results of 
which we can compare later with the above denominator and numerator zero 
previews. We give details for the case of eqn. (^) and quote our results for 
the other couplant equations. 

3 NUMERICAL SOLUTION OF THE [1|1] 

OPTIMIZED PADE QCD COUPLANT EQUA- 
TION OF MOTION 

Equation (|55D can be solved numerically for the couplant a, at any one chosen 
value of momentum Q and flavor number Nf, which means at any one chosen 
value of pi{Q, Nj) and p2{Nf). Then the value of Q and Nf can be changed 
and the numerical solution process repeated. We chose to keep Nf fixed at 
any one value of Nf in the range : < Nf < 16, while we varied Q over a 
wide range from very small to very large Q, and obtained a set of values of the 
(optimized) couplant a, that are solutions of eqn. at each chosen value 
of Q. In the process one can plot out a graph of these solutions of eqn. ([55|) 
against Q, for any one fixed value of flavor number Nf, to see how Pade QCD 
couplant, for a given flavor Nf, behaves over a wide range of Q values, from 
large Q ^ C)0, to small Q — ^ 0. Subsequently, the value of flavor number 
Nf can be changed and the entire process repeated to obtain a separate 



(a, Q) plot at the corresponding new value of Nf. Equation (|55|) was solved 
separately in this way for all integer flavor numbers Nf = 0, 1, 2, 3, 4, ... , 16. 
We point out immediately that our varying Q over the wide range Q —>■ oo,to 
Q ^ does not in any way imply that a Pade couplant solution necessarily 
exists in these momentum limits. As we explain further below, the procedure 
enables the system to pick its own cut-off momentum below which or above 
which a solution exists or does not exist. This will become fully clear below. 

We now give details of the exact method of numerical solution used, the 
method being the same for any one of the above Nf values. Our method was 
to regard the couplant a in eqn. (|55|) as a floating variable that is allowed to 
assume a wide range of values at any one fixed Q value (and also fixed Nf), 
such that for any one floating value of a, the left and the right hand sides of 
eqn. (|55| ) are in general not equal, meaning that the floating value of a, is in 
general not a solution of eqn. (|55|), and therefore not acceptable as a Pade 
QCD couplant. 

However, there exists some unique value (or multiple values) of the float- 
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ing variable a, (at a given Q and fixed Nf) for wliicli tlie left and tfie right 
liand sides of eqn. (^) are exactly equal, measured by a curve crossing a 
solution line. That unique floating value of a, becomes identifiable as our 
exact solution of eqn. (^). Thereafter, we change the value of Q and again 
keep this Q fixed at its new value, while the couplant a floats again over 
a wide range (and in very fine steps we chose as: Aa = 0.00001), until we 
again find that unique value of couplant a (crossing point) , for which the 
left and the right hand sides of eqn. (^) are exactly equal. In this way it 
was possible to obtain a set of couplant values a, that are exact solutions of 
eqn. ( |53| ) at given values of Q, all the time keeping to one fixed Nf value. 
For a different Nf value the entire search process is repeated, always starting 
from very low Q values, and moving up in suitable steps to higher Q values, 
or vice versa. The further details of the above manner of floating the cou- 
plant a, and constantly comparing when the left and the right hand sides of 
eqn. ( [55| ) are exactly equal, are what we describe next. 

First, we note that the double- valuedness of the optimized variable C2 in 
eqn ( pOD implies that eqn. (^) can be written out as two separate equations, 
each one of which can be separately solved numerically to yield its own 
(optimized) Fade couplant solution, together with a corresponding separate 
(a, Q) plot. The two separate equations we get out of eqn. (^) are: 



p,{Q,Nf)-P, = - + c\n 

Qi 



ka 

1 + cd — 2dPi 



(56) 



and 



PiiQ,Nf)-P2 = - + c\n 
a 



where: 



ka 



1 + cd- 2dP2 



(57) 



2 2 2c ^ ' 



and 



(5J 



with C2(+) and C2{—), being the two possible values of the optimized C2 
of eqn. (|UD. 

Next let us denote the left hand sides of eqns. (|56D and (^) by: 
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Xi = pi(Q,iV;)-Pi 



X2 = pi{Q,Nf)-P, 
Similarly, we denote the right hand sides of the same eqns. (|56|) and (|57|) 

by: 

ka 

1 + ca — 2aPi 
ka 

l + ca- 2aP2 

Then the two equations (^Up and (0) we want to solve numerically for 
a, can be written simply as 

Xi-X3 = Fi(a) = (60) 

X2 - X4 = F2(a) = (61) 

so that what we have to look for as the exact solutions of eqn. (|55D are 
the zeros of the two functions, Y'i(a) and 1^2 (a). We do this by assigning 
a wide range of values, in small incremental steps (Aa = 0.00001), to the 
floating couplant a, all at one fixed value of Q (and fixed Nf ). At any 
assigned value of a, we compute Yi{a) and 1^2(0). Then over the chosen 
wide range of a, (and at one fixed Q value and fixed Nj), we plot the two 
separate graphs of (Yi{a), a), and (^2(0), a)- The exact numerical solutions of 
eqns (^) and (|5^) we are looking for at a fixed Q, can now be read off these 
(Yi(a), a), and (1^2(0), «) plots, as the points on the couplant axis, where each 
curve crosses the couplant axis. Typical such plots and crossing points are 
shown in figs. |1] to ^ for the (Yi{a),a) solution, and figs. ^ and ^ for the 
(1^2(0), cl) solution. 

A simple MATLAB program we used for plotting out these (Yi(a),a), 
and (^2(0), a) curves, and reading off the crossing point values of a, on the 
couplant axis, is reproduced as Appendix A. The crossing point readings 
from the (Fi(a),a) plot, being solutions of eqn. ([56|), we have called the Yi 
component of Fade QCD. Similarly, we can call the crossing point readings 
from the (1^2(0-), a) plot, being solutions of eqn. (p7|), the Y2 component of 



X3 = - + cln 
a 



X4 = - + cln 
a 
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Figure 1: Sample crossing point profile on the Yl = line, Nf = 2 case 
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Graphical search for crossing point solution: Nf = 6 case 
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Figure 2: Sample crossing point profile on the Yl = line, Nf = 6 case 
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Graphical search for crossing point solution: Nf = 8 case 
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Figure 3: Sample triple point crossing profile on the Yl = line, Nf = 8 case 
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Graphical search for crossing point solution: Nf = 13 case 
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Figure 4: Sample crossing point profile on the Yl = line, Nf = 13 case 
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Graphical search for crossing point solution: Nf = 15 case 
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Figure 5: Sample crossing point profile on the Yl = line, Nf = 15 case 
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Graphical search for crossing point solutions: Nf=0 case 
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Figure 6: Triple point crossing configuration on the Yl = line, Nf = case 



22 



0.05 



0.04 



0.03 



0.02 



Graphical search for crossing point solutions: Nf = 12 case 
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Graphical search for crossing point solutions 
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Figure 8: Crossing point profiles on the Y2 = line 
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Figure 9: Crossing point profiles on the Y2 = line 
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Pade QCD. We find that the two solutions or components of Pade QCD are 
in general different and distinct from each other, implying immediately, that 
the [1|1] optimized Pade QCD couplant has at least two separate components 
(solutions), we shall from here onwards refer to as the Yi and Y2 component 
couplants of Pade QCD. 

4 FEATURES FOUND IN THE [1|1] OPTI- 
MIZED PADE COUPLANT EQUATION 

We now state the features we found with the above graphical and computa- 



tional analysis of the optimized [1|1] Pade QCD couplant eqn. (|55D. Later 
we point out which of these features persist as intrinsic features of the other 
Pade QCDs we also analysed by the same graphical method. 



4.1 The Multiplicity Structure in the Yi Pade Compo- 
nent Solution 

One striking feature we found in the course of studying the above (Yi{a), a), 
and (1^2(0), a) Pade QCD plots, is that while the {Y2{a), a) plot is simple and 
structureless for all flavor numbers Nf, eqn. ( |56[ ) and its (^1(0), a) plots have 
in general three different but simultaneous crossing points on the couplant 
axis, for a given Q value, implying that eqn. ( |5B| ) has in general three different 
solutions for the same Q value and the same flavor number, and this feature 
holds for each flavor number : < Nf < 16. Typical such (simultaneous) 
triple crossing points can be seen in figs. to |^. The feature exists for every 
flavor number < Nf < 16. 

The first graphical indication of the existence of this triple multiplicity 
solution of eqn. (|56|) is the observation that as Q increases from below, the 
entire profile of the {Yi{a),a) plot rises upwards, such that beginning with 
the plot (at low Q) cutting the couplant axis at only one distant point to the 
far right hand side of the couplant axis, the curve later rises enough (at higher 
Q values) to begin to cut and cross the couplant axis simultaneously at three 
distinct points. We found that this begins to occur when Q has attained a 
certain minimum or threshold value we have denoted by Qmin, and that this 
minimum Q value exists as a sharply defined threshold point for each flavor 
number Nf . Significantly however, the actual value of Qmin, (see Table ^, 
differed from one flavor number system to another. These rising profiles of 
the (Fi(a),a) plots and their triple crossings for Q > Qmin are shown in 
figs. [1^ to |l^ for some flavor numbers, the pattern being however the same 



26 



for all flavor numbers. Figures to |1^ show sample profile positions when 
one is at exactly the threshold (critical) point Q = Qmin for these rising 
profiles. 

We observed also from the same {Yi{a),a) profile plots that the triple 
crossing point feature does not persist for all Q > Qmin- Rather, as Q rises 
higher and higher above Qmin, the triple crossing suddenly ceases and we 
noted that this occurs when Q has attained some upper cut-off value we 
denoted hj Q = Qmax- This upper cut off value of Q was found to be as 
precise and sharply defined for each Nj flavor system as the lower threshold 
point Q = Qmin- But, again as in the Qmin case, the value of Qmax, (see 
Table differed from one Nf flavor system to another. The upper parts 
of Figures to show sample profiles and configurations at the cessation 
point of triple crossings when Q = Qmax- 

By way of comparison, we show in figs. EUI to 521 the profiles of the 
(12(0); 0,) plots for some flavor numbers. None of these has a triple crossing 
structure, but only one solution or crossing point for eqn. (^) at any one 
value of Q, and for all flavor numbers. 

Focusing now on these crossing point solutions of eqns. (|56D and (^7]), 
we can denote by 01,02,03, the general triple point crossing solutions of 
eqn. (|56D, and by 04 the one solution (crossing point) of eqn. (|57D. Here 
Oi is the extreme left crossing point in the (Yi(o),a) plot, 02, is the middle 
crossing point, while 03, is the extreme right crossing point, all in the same 
(Y'i(o),o) plot of any given (fixed) Nf system. Together the four solutions 
01,02,03,04 imply that the Fade approximant QCD couplant o of eqn. (|55|) 
or eqn. (|ID, has a multiplicity structure of four distinct component couplants 
or color force solutions, we can read off as crossing points on our Yl = 
and Y2 = solution lines. This multiplicity of solutions of the Fade 
couplant equation (|55D, exhibited in the momentum band Qmin < <5 < Qmax, 
represents our first explicit finding concerning the features of Fade QCD 
couplant analyzed by our above computational and graphical method. 



4.2 Features of the ai{Yi) , a2{Yi) , a^{Yi) , and 04(12) Fade 
couplant Solutions. 

We investigated next the behavior of each Fade component couplant solution, 
by plotting its graphical Fade crossing point values against the momentum 
Q of the profile plot at which the crossing point value was read off. We found 
the following features. 
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Figure 10: Yl rising profiles, Nf = 2 case 
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Figure 11: Yl rising profiles, Nf = 9 case 
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Figure 13: Yl rising profiles, Nf = 15 case 
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Graphical search for crossing point solutions: Nf = 16 case 
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Profile at onset.and at cessation of triple point crossing: Nf = case 
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= case 
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Figure 16: Critical momentum profile configurations on the Yl = line, Nf 
= 2 case 
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Figure 17: Critical momentum profile configurations on the Yl = line, Nf 
= 8 case 
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Profile at onset, and at cessation of triple point crossing: Nf = 10 case 
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Figure 18: Critical momentum profile configurations on the Yl = line, Nf 
= 10 case 
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Profile at onset,and at cessation of triple point crossing: Nf = 15 case 
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Figure 19: Critical momentum profile configurations on the Yl = line, Nf 
= 15 case 
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Figure 20: Y2 rising profiles, Nf = case 
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Y2 rising profiles: r^=l4o^^GeV 
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Figure 21: Y2 rising profiles, Nf = 4 case 
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Figure 22: Y2 rising profiles, Nf = 13 case 
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4.2.1 Features of the Fade component couplant solution 

From fig. |I| we see tliat tlie as is tlie only solution or crossing point of eqn. (|56|) 
in the low energy momentum region Q < Qmin, and for any given flavor 
number Nf, particularly in the range, < Nf < 8. For 9 < Nf < 16, the 



indication from flgs. |TT] to [Tj as well as figs. ^ , ^ and is that the 03 
is already infinitely large even before we reach the point Q = Qmin- If we 
combine this with the features of the 03 shown in fig. for the < Nf < 8 
cases, the indication is that while the 03 does exist in general in the region 
Q < Qmiri, it has the form of a Landau type pole singularity for a Fade 
beta function or couplant. This is re-inforced by the fact seen in fig. ^ that 
even for the lowest A'^^ = case, the 03 already rises sharply towards very 
large values at Q = 0.705 GeV, in the process still cutting off access into the 
infra-red region Q < 0.705 Gev. 

The 03 exists in the medium energy region Qmin < <5 < Qmax but with 
a value that decreases progressively with increasing momentumQ. It finally 
cuts off at the characteristic higher momentum Q = Qmax such that for 
all Q > Qmax, CLs = 0. Because this second Fade couplant bifurcation 
point occurs well inside the normal FQCD region of QCD, Q > Aqcd, with 
Qraax ^ ^QCD, the possibility arises that if the above Fade couplant structure 
represents physical QCD reality, the Fade (03,02) structure can affect sub- 
stantially even relatively high energy QCD processes, by way of distinctively 
non-perturbative contributions. 

4.2.2 Features of the ai Fade component couplant solution 

Examining the ai Fade component couplant solution in the same way, we 
observe from fig. ^ as well as the upper halfs of figs. |TB| to |r^, that the ai 



solution (crossing point) is the only Fade QCD present in the very high 
energy region Q > Qmax where it also has very small values for any fiavor 



number, such that for Q ^ 00, Oi — > 0, shown in fig. 

The ai couplant solution exists in the intermediate energy region: Qmin < 
Q < Qmax attaining its highest but still moderate value at Q = Qmin from 
where it decreases progressively towards zero as Q ^ 00. The oi solution 
does not exist in the low energy region < Q < Qmin- 

The indication from these features is that the Fade oi component solution 
can be identified as the asymptotically free purely FQCD component color 
force of QCD we started with in eqns. (P and (^. We shall so identify the Oi 
component Fade solution. Then the fact that the ai solution does not exist 
in the region < Q < Qmm is seen as a correct refiection of the denominator 
zero singularity of the Fade beta function eqn. (|20|) and its attendant Kogan- 
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Shifman 0] type behavior, with our Qmin directly identifiable with the /ic cut 
off momentum of Kogan and Shifman, and of Elias et. al. [|l], 0|. 

Based on these identifications, the question of whether infra-red scenario 
I (frozen couplant), or scenario II (bifurcated infra-red attractor point) holds 
at any fiavor in Pade QCD and perhaps in real QCD, boils down in our 
approach to the question of the extent the momentum gap < Q < Qmm 
remains finite and free of our (01,02) bifurcating components, in various 
flavor states of QCD. We will find that our method allows us to answer this 
question explicitly later. 



4.2.3 The 02 Pade component couplant solution 

Plotting the 02 crossing point values against momentum Q, we find the be- 



havior shown in fig. |2^. That is, the 02 component color force solution is 
found to exist only in the intermediate energy region Qmin < Q < <5max for 
any fiavor number Nf, but vanishes or cuts off sharply at the two critical 
upper and lower momentum points Q = Qmax and Q = Qmin- 



From fig ^ we see that inside its intermediate energy domain of operation 
and existence, the 02 couplant rises in value from a lowest value at Q = Qmin 
to a highest value at Q = Qmax- This manner of behavior suggests we identify 
our Pade 02 component with the upper branch of the bifurcated structure 
found by Elias et. al. 0] and by Kogan and Shifman This is reinforced 
by the observation from our profile plots shown variously in figs. |l] to |19] that 
the 02 crossing point solution merges with and coincides exactly in value with 
the ai couplant solution at the point Q = Qmin which can be regarded as a 
bifurcation point of the ai into 02. Then the infra-red region < Q < Qmm, 
is totally free of the two bifurcated component couplants, ai and 02, again 
agreeing fully with the behavior found by Elias et. al. fl], 0, and by Kogan- 
Shifman 0. 

However, unlike the Elias et. al. structure which did not determine the 
progression of this 02 upper branch, our graphical crossing point solution 
method shows the 02 clearly as terminating and merging exactly with the 
03 crossing point solution at the point Q = Qmax, where 02 = 03, and both 
couplants thereafter disappear together for all Q > Qmax- One may say that 
the 02 bifurcates further into the 03 which then runs towards the infra-red 
region but soon rises sharply to very large values in a manner suggestive 
of a Landau pole behavior. This 03 presence was not discernible from the 
numerator and denominator zero analysis of Elias et. al. 0]. 



43 




Figure 24: Pattern of variation with momentum and flavor, of isolated Pade 
Qi component couplant solution 
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4.2.4 The 04(1^2) Pade component couplant solution 

The 04(12) crossing point solution, arose from the separate eqn. (|57|). As such 
it emerges for now as a lone star Pade component color force, not connected 
to the 01,02,03 components of eqn. Its features are shown in fig. 

where it is seen that it has no structure or bifurcation but behaves smoothly, 
falling asymptotically in value from a large value (04 00) near Q — > 0, 
to the value (04 — 0) as {Q — > 00). The exact meaning of the 04 couplant 
solution is not yet clear to us, and is being investigated. 

4.3 Spiral Chain-like Structure found in Pade QCD 

Arising from the fact stated above that while the 01,02,03 are independent 
solutions of Pade eqn. (|56|), they are however joined together at the two 
critical (bifurcation) momentum points Qmin and Qmax? in the manner given 
by: 

01 = O2 7^ 0, at : Q = Qmin, 

and 

02 = O3 7^ 0, at : Q = Qmax- 

we find that a combined plot of the individual variations (oi, Q); (02, Q); (03, Q) 
gives one unbroken chain-like spiral structure, shown in figs. to ^ for var- 
ious flavors. The 04 component has also been plotted in. 

These chain-like spiral structures together with our computed values of 
the bifurcation point critical momenta Qmin and Qmaxc shown in Table |^, 
lead us to explicitly answer the question, for which flavor states the infra-red 
scenario I or scenario II holds. Additionally, we have also shown in Table |^, 
the values of the Pade couplants oi((5min) and 03(Qmax) at these critical 
cut-off momentum points. They all enable us to answer the same question 
definitively. Before discussing the question however, we first show that the 
above features of the optimized [1|1] Pade QCD are also shared to a large 
extent by the other Pade QCDs, [1|1], [2|1] and [1|2] that we also analyzed. 

5 ANALYSIS AND FEATURES OF THE NON- 
OPTIMIZED [1|1], [2|1] AND [1|2] PADE QCDs 

. Using exactly the same computational and graphical procedures described 
above for the optimized [1|1] Pade QCD, we also analyzed the non-optimized 
[1|1], [2|1], and [1|2] Pade QCD couplant equations (g5|) , (^ and (|31|). We 
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Variation of a4 witin Energy and Flavor 
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Figure 26: Pattern of variation with momentum and flavor of Pade 04 com- 
ponent couplant solution 
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Figure 27: A plot showing a spiral multi-component structure of Pade QCD: 
Nf = 3 case 
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QCD Spiral structure: Nf = 3 




200 400 600 

Energy Q (Gev) 

QCD Spiral structure: Nf = 3 



800 











I 










a3 




















Qmax 


\a1 
























.Qmin 

















200 400 600 
Energy Q (Gev) 



800 



0.0824 



U3 

-§.0.0823 

o 
O 

I 0.0822 

c 

o 

Q. 

E 

O 0.0821 



QCD Spiral structure: Nf = 3 



0.082 





0.1 








CO 





09 


Q. 






13 

o 





08 


O 











07 


0) 






c 
o 





06 


Q. 






E 
o 





05 ■ 


O 











04 - 


> 











03 




















a3 























200 400 600 

Energy Q (Gev) 

QCD Spiral structure: Nf = 3 



800 











{ ^ 

















































200 400 600 
Energy Q (Gev) 



800 



Figure 28: A plot showing a multi-component couplant structure in Fade 
QCD: Nf = 3 case 
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Spiral Multi-component structure of QCD: Nf = 3 
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Figure 29: A plot showing a spiral multi-component structure in Pade QCD: 
Nf = 3 case 
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Spiral Multi-component structure of QCD: Nf = case 







1 

Nf = case 






ZOOM used, axis([-100 1 .5e+3 0.01 0.12]) 


















a3 




r 


a2 




a2 


Qmax 










\ a1 














a1 


1 


Qmin 


a4 

1 



500 



1000 



1500 



Energy Q (GeV) 



Figure 30: A plot showing a spiral multi-component structure in Pade QCD: 
Nf = case 
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Spiral Multi-component structure of QCD: Nf = 1 case 
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Figure 31: A plot showing a spiral multi-component structure in Pade QCD: 
Nf = 1 case 
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Figure 32: A plot showing a spiral multi-component structure in Pade QCD: 
Nf = 2 case 
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QCD Spiral structure: Nf = 2 
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Figure 33: A plot showing a multi-component couplant structure in Fade 
QCD: Nf = 2 case 
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Figure 34: A plot showing a spiral multi-component structure in Pade QCD: 
Nf = 4 case 
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Spiral Multi-component structure of QCD: Nf = 6 case 
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Figure 35: A plot showing a spiral multi-component structure in Pade QCD: 
Nf = 6 case 
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Spiral Multi-component structure of QCD: Nf = 8 case 
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Figure 36: A plot showing a spiral multi-component structure in Pade QCD: 
Nf = 8 case 
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QCD Spiral structure: Nf = 8 
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Figure 37: A plot showing a multi-component couplant structure of Fade 
QCD: Nf = 8 case 
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Figure 38: A plot showing a spiral multi-component structure in Pade QCD: 
Nf = 10 case 
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Figure 39: A plot showing a multi-component couplant structure in Fade 
QCD: Nf = 10 case 
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Figure 40: A plot showing a spiral multi-component structure in Fade QCD: 
Nf = 15 case 
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QCD Spiral structure: Nf = 15 
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Figure 41: A plot showing a multi-component couplant structure of Pade 
QCD: Nf = 15 case 
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Figure 42: A plot showing a spiral multi-component structure in Pade QCD: 
Nf = 12 case 
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QCD Spiral structure: Nf = 12 
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Figure 43: A plot showing a multi-component couplant structure of Fade 
QCD: Nf = 12 case 
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found essentially the same Pade QCD features, with small differences which 
we summarize as follows. 

1. The 04 component solution found in optimized Pade QCD is not exhib- 
ited by any of the three non-optimized Pade QCDs. As will be recalled, 
the 04 solution arose from the quadratic nature of the optimized C2 co- 
efficient given in eqns. (^Of ) and dSOf ). This C2 quadratic feature is not 
present in the non-optimized [1|1], [2|1], [1|2] Pade QCDs, and as such 
the 04 solution does not arise in these cases. We can therefore for the 
moment play down the existence and role of the Pade 04 component 
solution. 



2. However, the triple crossing point feature and hence the ai, 02, 03 com- 
ponent structure, are present in exactly the same form in all the four 
Pade QCDs we analyzed, optimized or not. Correspondingly, the spiral 
chain-like structure exists in all of them and can be taken as intrinsic 
characteristic feature of a Pade QCD. 

3. We found that while each Pade QCD at a given flavor, is characterized 
by two critical or cut off momenta Qmin and Qmax, the actual values 
of these momentum pairs differed substantially from one Pade QCD 
to other. The values found for example, for the non-optimized [1|1] 
Pade QCD are shown in Table ^ to be compared with Table ^ for the 
optimized [1|1] Pade QCD. The [2|1] and [1|2] paired critical momenta 
differ just as much. The trend of the variation is not yet clear but is 
being separately studied. 

4. In contrast, the values of the critical (cut-off) couplants ai{Qmin) and 
c^slQmax) exhibit some consistency, especially in high flavor states seen 
in Tables ^ and ^. 

5. We found that the optimized [1|1] Pade QCD has the feature of resolv- 
ing clearly, the spiral nature of the Nf = 6, 7, 8 flavor states shown in 
flgs. ^ and compared to the non-optimized [1|1] Pade QCD where 
the couplant structure of the same Nf = 6, 7, 8 flavor states, is not so 
resolved, as can be seen in flg. 



The above small differences apart, we affirm that our Pade QCD of various 
orders, optimized or not, has the characteristic feature of triple multiplicity 
of solutions (01,02.03); two characteristic bifurcation points, Qmin and Qmax; 
and a spiral chain-like structure continuously connecting the three component 
solutions, (oi, 02,03). Using flgs. ^ to ^ of the [1|1] NNLO optimized Pade 
QCD as typical of these behaviors and properties of a Pade QCD, we now 
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consider the question of which of infra-red scenario I or II holds at a given 
flavor. 



6 THE QUESTION OF IR SCENARIO I OR 
SCENARIO II IN FADE QCD 

As already stated, in our approach, the question of whether we have infra-red 
scenario I (IR stable fixed point and frozen couplant) or IR scenario II (IR 
attractor point and bifurcated couplant) is directly decided by whether the 
momentum gap < Q < Qmin exists and is finite. Where this is the case, 
the IR scenario II necessarily holds, ruling out scenario I. Conversely, if the 
momentum gap does not exist for a given flavor, the IR structure of that 
flavor state is necessarily of scenario I type, and scenario II becomes ruled 
out. We now examine our spiral plots as well as our Table § to see which 
scenario holds. We find as follows: 

1. The case < iV/ < 5 

From our spiral plots figs, to it is clear that IR scenario II is what 



operates and not scenario I. The same conclusion is drawn from our 
Tables |^ and ^, where it is seen that the momentum gap < Q < Qmin 
is far from being zero. These findings agree with those of Elias et. al. 
who found that scenario II holds for all < Nf < 5, regardless of which 
Fade approximants, [1|2], [2|1], [1|3], [3|1], or [2|2] they used. 

2. The case Nf = 6, 7, 8 



From our spiral plots shown in figs. |35| and ^ it is also clear that 
IR scenario II holds for Nf = 6, 7, 8, and not scenario I. The Qmin 
values shown in Tables |^ confirm this further. The momentum gap 

< Q < Qmin is far from being zero in all these cases. 

3. The case 9 < iV/ < 16 

For these flavor states, Elias et. al. by their method, found scenario 

1 IR behavior. In contrast, our method shows the opposite scenario II 
behavior for all 9 < Nf < 16 as for all < Nf < 8. Our finding of 
scenario II for all 9 < Nf < 16 is seen clearly in figs. |3^ to ^ for the 
cases of Nf = 10, 12, 15, chosen for illustration. All the other flavor 
states Nf = 9, 11, 13, 14.16, have exactly the same scenario II behavior. 

We can explain this particular difference in the two results, and estab- 



lish that our finding is the correct Fade position. Looking at figs. 38 



to ^ as well as Tables ^ and |^, it is seen that the infra-red attractor 
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points exist at progressively lower and lower momentum in these high 
flavor states, and only a method designed to follow this trend can es- 
tablish the continued existence of the scenario II behavior. The Elias 
et. al. method would appear not equipped for this. 

Even with our method, looking only at our figs. PSI, EDI, and ff^, one 



would gain the impression that our ai and 02 Fade component cou- 
plants which are the exact analogue of Elias et. al. bifurcated cou- 
plants at /i = fic, stretch on to zero momentum point, Q = /i = in 
which case scenario I could be concluded. 

However, if one follows these two branch Fade component couplants 
Oi and a2, to sufficiently low momentum regions, shown typically in 
our figs. ^ and ^ one sees unambiguously, that these two component 
couplants still meet and turn back at some finite critical momentum 
Qmin = yUc 7^ 0, leaving a clearly visible gap, < Q < Qmin, where both 
the ai and 02 are totally absent. This gap is present and completely 
visible in all cases 9 < Nj < 16, leaving no doubt that scenario II is 
intrinsically what holds in these Fade fiavor cases as for the < Nf < 8 
cases. 

It just happens that for the higher fiavor states, the spiral (bifurca- 
tion) point where ai and 02 meet, shifts to lower and lower momentum 
point, and unless this shifting is followed from fiavor to fiavor, one 
can conclude erroneously that scenario I operates in these higher fiavor 
states. 

Against the above, we can affirm that our results are actually consistent 
with those of Elias et. al. and that the features we found and analyzed 
in our Fade QCDs, are the intrinsic features of Fade QCD of any order 
[1|1], [2|1], [1|2], [3|1], [1|3] or [2|2]. However our graphical method of analyz- 
ing these Fade QCD features, has the advantage of showing clearly that the 
infra-red scenario II behavior is the intrinsic behavior of Fade QCD in all 
fiavor states < Nf < 16, and not IR scenario I. 

We are led also to firmly identify our critical momentum point Qmm = f^c 
as fundamentally an infra-red attractor point or a pole singularity of Fade 
QCD beta function, for all fiavor states, < Nf < 16. This means that the 
critical Fade QCD couplants ai(Qmin) shown in Table || are fundamentally 
infra-red attractor points and not IR stable fixed points, although for suffi- 
ciently high fiavor Nf, the attractor behavior (point) can be mistaken for a 
IR stable fixed point or frozen couplant behavior. 
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7 COMPARISON WITH TRUNCATED PQCD 
INFRA-RED FIXED POINT RESULTS 



Having settled the above question, we can check for further consistency of 
our Pade QCD compuatations and findings by comparing them with the 
earher work of Banks and Zaks 0; and of Stevenson et. ah 0, 0, who 
used purely truncated PQCD beta function to search for infra-red (IR) fixed 
points in PQCD. Banks and Zaks |^ as well as Stevenson |^ analyzed the 
NLO truncated PQCD of eqn. (H) and predicted PQCD IR fixed points shown 
against flavor numbers in our Table ^ Separatelyly, Stevenson et. al. 0, 
^, |], using optimized NNLO truncated PQCD of eqn. investigated the 
same IR fixed point structure of PQCD. The IR fixed points they found 
are shown in the same Table |[ We can now compare these results with 
our fli values computed from our [1|1] optimized (NNLO) Pade QCD. These 
ai(Qmin) values are also shown side by side in Table |[ We now observe as 
follows. 

1. For low flavors Nf, there is no agreememt between our ai values and 
the NLO or NNLO IR fixed point values of the above authors. This is 
consistent with the earlier stated fact that the point (Qmin, ai(Qmin)) is 
not an IR stable fixed point of QCD but a pole singularity (IR attractor 
point) of the beta function. 

2. For very high Nf = 16, 15, 14, 13.., our values of ai agree almost exactly 
with the NLO and NNLO IR stable fixed point values of the above 
authors. This agreement is evidence that our Pade QCD computations 
and graphical solution method are correct and reliable. However, the 
agreement comes not because our oi {Qmm) point has suddenly become 
IR stable fixed point of scenario I, but from the spiral structure shifting 
to lower and lower momentum points discussed earlier. As a result of 
this shifting, the IR attractor point ai(Qmin) becomes mimiced more 
and more effectively by an IR stable fixed point of a seemingly frozen 
PQCD (or Oi) couplant. But when one follows this bifurcation point 



shifting to sufficiently low energy as shown typically in our figs. ^ 



and ^ one soon sees that ai and 02 are still separate or bifurcated 
component couplants and that the PQCD couplant ai did not freeze 
at all. 

We conclude finally that our Pade QCD graphical computations and 
findings, can be relied upon, being consistent with the independent Pade 
beta function numerator /denominator zero analysis of Elias et. al.; and also 
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consistent with the direct PQCD IR computations of Banks and Zaks, and 
Stevenson et. al. 

8 SUMMARY AND CONCLUSIONS 

We have presented a rehable and consistent graphical computational method 
of analyzing the infra-red structure of Fade approximant QCDs. We have 
shown with the method that the intrinsic behavior of Fade QCD is one 
in the FQCD bifurcates at some critical Kogan-Shifman type momentum 
f^c = Qmin, leaving the infra-red region < fj^ < fj^c totally decoupled from 
the FQCD dynamics. We found however that besides this FQCD bifurcation, 
there exists a second bifurcation by the upper branch 02 couplant, into a third 
branch or component couplant 03, at another critical momentum Qmax where 
in all cases, Qmax > <5min- The three components or branches (ai, 02, 03) are 
arranged in a chain-like spiral structure. The 03 in general runs into the 
original infra-red region ^ < but soon rises sharply to very large values in 
a manner suggestive of a Landau pole behavior of the Fade couplant at its 
second bifurcation point. The extent to which the above Fade QCD features 
and structure can be taken over as those of real physical QCD is a matter 
for further studies which we are now following up. 
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A (Appendix A): The exact MATLAB Pro- 
gram (QCDS.m) used by the author to 
solve numerically, the optimized SMQCD 
eqn. ( [55l ) of this paper 

% The footnotes given below are to guide the actual execution of the pro- 
gram; they are not part of the program. 
% The program now begins. 

clear all; 

= 2; % see footnote No. 1 below. 
Rh2 = -10.91120013471925; % see footnote No. 2 below. 
Q = 1.0; % see footnote No. 3 below. 

L = 0.230; % L = Aqcd = 0.230 GeV, kept constant all through the SMQCD 
computations. 

rl = 1.9857 — 0.1153 * A^; % rl is a constant all through the (optimized) 
SMQCD computations. 
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b = (33 — 2 * N)/6; % 6 is a constant all through the SMQCD computations, 
c = (153 - 19 * iV)/(66 - 4 * A^); % c is a constant all through the SMQCD 
computations. 

Rhl = b* log{Q/L) + c * log{2 * c/b) — rl; % see footnote No. 4 below. 

D = 8*c*c — 4* Rh2] % D varies as we vary input Rh2. 

PI = 3 * c/2 + sqrt{D) /2; % PI varies as we vary input Rh2. 

P2 = 3* c/2 — sqrt{D) /2; % P2 varies as we vary input Rh2. 

XI = Rhl - PI; % XI varies as we vary Q and Rh2. 

X2 = Rhl - P2; % X2 varies as we vary Q and Rh2. 

a = 0.02 : 0.00001 : 0.12; % see footnote No. 5 below. 

Fl = 1.0 + c * a — 2 * a * PI; % Fl is just a notation. 

F2 = 1.0 + c * a - 2 * a * P2; % F2 is just a notation. 

HI = ((c * a)./Fl); % note the space and the dot before /Fl 

H2 = ((c * a)./F2); % note the space and the dot before /F2 

X3 = (l./a) + c * log{abs{Hl)); % note the space and the dot before /a 

X4 = (l./a) + c * log{abs{H2)); % note the space and the dot before /a 

Yl = XI - X3; 

Y2 = X2 - X4; 

plot(a, Yl); grid on % see footnote No. 6 below. 
xlabel(' Optimized couplant') 

ylabel('Yl Deviation from optimization solution (crossing point)') 
title('Graphical search for crossing point solution of the Optimization Eqn. 
for N =2') 

% use next a suitable axis command a.xis{[xiX2yiy2]) , to zoom closely into a 
crossing region to be read off. % This zooming in with the axis command 
needs to be used all the time, and with great skill, to obtain accurate read- 
ings. % The MATLAB program ends here. 

% FOOTNOTES begin (not part of the program). 

% Footnote No. 1: After varying Q over all desired ranges at one fixed 
flavor number N, we change N by hand, by simply inserting any other value: 
N = 0, 1, 2, 3, 16. 

% Footnote No. 2: The Rh2 values we used were those computed pre- 
viously by Mattingly and Stevenson ^ Explicitly, these p2 values 
are: 

% Nf = 0,p2 = -8.410032589173554; Nf = I, P2 = -9.996607149709793; 
% Nf = 2,p2 = -10.91120013471925; Nf = 3, p2 = -12.20710268197531; 
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% Nf = 4,p2 = -13.90995802777778; Nf = 3,p2 = -15.49181836878120: 



%Nf = 


6, P2 = - 


-17.66469557734694; Nf = 


= 7, P2 = 


-19.78668878025239; 


%Nf = 


8,P2 = - 


-22.74511792421761; Nf = 


= 9, P2 = 


-25.96983971428571; 


%Nf^ 


10,P2 = 


-30.64825148592373; Nf 


= 11,P2 


= -36.70527878387512 


%Nf = 


12,P2 = 


-46.58505774333333; Nf 


= 13, p2 


= -63.56012999671129 


%Nf^ 


14, P2 = 


-101.9145678055556; Nf 


= 15, P2 


= -229.8874551066667 


% Nf = 


16, P2 = 


-1724.404563921111. 






% Footnote No. 


3: Start in general with Q = 1.0 GeV (for any flavor). 



and after observing the SMQCD profile pattern, you vary Q upwards or 
downwards in whatever steps (large or small) you desire, but state the value 
of Q in GeV always, e.g. Q = l.Oe - 305 for Q = lQ-^°^ GeV. 

%Footnote No. 4: For 9 < N < 16, you must replace c*log{2*c/b) — rl; 
by c * log{abs{2 * c/b) — rl; This is because c changes sign as we enter the 

> 9 or (c~) phase of QCD, so we need to use the absolute value. This 
apart, there is no other difference with the {c'^) or < < 8 computations. 
Notice however, that the SMQCD profile now turns upside down. 

% Footnote No. 5: The two end points: a — 0.02 and a = 0.12 are in 

general not fixed for any one N value but can be chosen differently or varied 
in the course of any one computation. The optimal choice to suit an already 
chosen Q value will be found by experience. In general, when Q is very 
low, we need to make the terminal point of the a range very large, in order 
to be able to read at all the as{Yi) crossing point. However, whatever the 
range chosen for a, the same incremental step Aa = 0.00001 of the floating 
couplant should be used throughout these computations to achieve the same 
degree of accuracy of crossing point readings (solutions). Also for any one Q 
value the chosen range of a should not be too large, otherwise the computer 
goes into a time consuming loop and floods its memory unduly. 

% Footnote No. 6: After plotting the (a, Yl) curves for all desired Q val- 
ues, we replace the command: plot (a, Yl), with the command: plot (a, Y2). 
Then we vary Q again over all desired values (with N still kept flxed) , in order 
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to plot out the (a, Y2) profiles and read off the crossing points (solutions). 

% Footnote No. 7: Corresponding changes are to be made by hand in 
the ylabel and title commands. 
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Table 1: Values of the denominator zeros (a^), and the numerator zeros (a„) 
of the Pade beta function eqn. (|2^) computed as a function of flavor number 



r lavoi iNuniDei i\ t 


uenomniatoi zeio 


iNumeiaLOi zeio 


U 


U.ZoOD 


U.o40O 


i 


U.ouy4 


U.yooy 


n 
L 


n Q /I QQ 

U.o4oo 


1 (V7 A n 

i.U r 4y 


q 
O 


U.oy /o 








2 2780 


5 


0.8549 


-10.9638 


6 


- 3.2000 


-0.8058 


7 


- 0.2280 


-0.2036 


8 


-0.0063 


- 0.0063 


9 


0.0.0799 


0.0763 


10 


0.1288 


0.1089 


11 


0.1625 


0.1150 


12 


0.1886 


0.1056 


13 


0.2108 


0.0873 


14 


0.2307 


0.0640 


15 


0.2494 


0.0384 


16 


0.2676 


0.0126 
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Table 2: Values of Qmin and Qmax in optimized [1|1] Padc QCD, tabulated 
here as a function of flavor number Nf , together with the corresponding 
values of ai((5mm) and a3((5max) at these spiral (bifurcation) points. 



Flavor Number Nf 


Qmin (GcV) 


'^1 (Qmin) 


Qmax (GcV) 







25.725856708 


0.063660 


1014.5456575 


0.074680 


1 


30.04865615 


0.065450 


722.2234986368 


0.076230 


2 


33.44905052 


0.068430 


907.79549390 


0.079170 


3 


37.90878530 


0.071720 


679.6977179750 


0.082210 


4 


43.8476891550 


0.075430 


1763.271161880 


0.085340 


5 


49.2803656516 


0.080480 


435.1917430243 


0.089580 


6 


56.6158149 


0.086260 


316.669113500 


0.093780 


7 


60.714770460 


0.094350 


184.198105110 


0.099280 


8 


56.96621400 


0.103870 


59.14194290 


0.104190 


9 


5.1612395010 


0.10950 


47.126555910 


0.117200 


10 


6.858398820(10-2) 


0.111920 


73.2831370350 


0.133120 


11 


3.742331065(10-=^) 


0.111070 


276.8901326 


0.154840 


12 


1.3532233520(10-6) 


0.103080 


1.019000784(10+^) 


0.180700 


13 


8.5876890(10-15) 


0.087960 


1.771666446(10+^) 


0.214860 


14 


3.540839347(10-^9) 


0.065770 


1.10405153(10+22) 


0.256120 


15 


6.26397660(10-1°^) 


0.039650 


2.17827906(10+89) 


0.310460 


16 


< 10-307 


0.01280 


> 10+308 


0.380525 
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Table 3: Values of Qmm and Qmax, in the non-optimized [1|1] Fade QCD. 
The gap between Qmin and Qmax closed up so fast after Nf = 5, that we 
could not resolve or read off the crossing point values for Nf = 6,7 and 8 in 
this non-optimized case. 



Flavor Number 
Nf 


(GeV) 


ai(Qmin) 


Qmax 

(GeV) 


a3(Qmax) 



1 
2 
3 
4 
5 
6 
7 


5.6468615736(10+^) 
4.4044582296(10+6) 

3.2521338660(10+5) 
2.1002809749(10+^) 
1.0134866543(10+''^) 
2.53433811590(10+^) 


0.009440 
0.011590 

0.014690 
0.019640 
0.02911 
0.058180 


2.2356088581(10+^) 
2.6026228349(10+^) 

1.40497854250(10+6) 
1.0761648724(10+5) 
8.5575884074(10+=^) 
2.351206362(10+2) 


0.009650 
0.011890 

0.015140 
0.020350 
0.03047 
0.062790 


8 
9 

10 
11 
12 
13 
14 
15 
16 


7.928886017(10+12) 

2.083364761(10+5) 
2.064339838(10-1) 

5.940796088(10"^) 
3.1465791130(10-22) 

2.617118166(10-56) 
4.098229275(10-168) 

< 1.0(10--^"') 


0.012690 
0.026420 
0.043050 

0.062140 
0.075920 
0.069880 
0.043470 

0.013230 


1.5839622507(10+13) 

3.5918756304(10+6) 

5.9926835830(10+3) 

6.6365417260(10+1) 
8.845270086(10-1) 
3.096682427(10-3) 
2.848034499(10-8) 

1.101656210675(10-3«) 


0.012790 
0.027450 
0.048350 

0.083840 
0.154850 
0.3321469 
0.99750 

9.631880 
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Table 4: Comparative values of the IR fixed points of PQCD from NLO 
& NNLO perturbative QCD calculations, and the values of ai((5min) and 
03(Qmax) shown earlier in Table ||, computed from the optimized [1|1] Fade . 



Flavor Number 


Qmin 


Ctl(Qmin) 


a*,^iPQCD) 




Nf 


(GeV) 


op. [1|1] 


at NLO 


at NNLO (optimized) 





25.725856708 


0.063660 


0.4112 


0.313284 


1 


30.04865615 


0.065450 


0.3863 


0.280270 


2 


33.44905052 


0.068430 


0.3614 


0.2634796 


3 


37.90878530 


0.071720 


0.3364 


0.244217 


4 


43.8476891550 


0.075430 


0.3115 


0.224065 


5 


49.2803656516 


0.080480 


0.2866 


0.208085 


6 


56.6158149 


0.086260 


0.2617 


0.190693 


7 


60.714770460 


0.094350 


0.2368 


0.176191 


8 


56.96621400 


0.103870 


0.2118 


0.160361 


9 


5.1612395010 


0.10950 


0.1869 


0.146072 


10 


6.858398820 x lO'^ 


0.111920 


0.1620 


0.1304388 


11 


3.742331065 x 10"=^ 


0.11107 


0.1371 


0.1150355 


12 


1.3532233520 x 10"^ 


0.103080 


0.1121 


0.0979828 


13 


8.5876890 x lO^^^ 


0.087960 


0.0872 


0.0797984 


14 


3.540839347 x 10"^^ 


0.065770 


0.0623 


0.05940013 


15 


6.26397660 x 10^^°^ 


0.039650 


0.0374 


0.03688832 


16 


< 10-307 


0.0128 


0.0125 


0.01248992 
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